data2<-read.table("smoothtest_merge2.txt",sep="\t",head=T)
merge2<-read.table("smoothtest_merge2.txt",sep="\t",head=T)
merge4<-read.table("smoothtest_merge4.txt",sep="\t",head=T)
dim(merge4)
dim(merge2)
rownames(merge2)<-as.character(merge2[,1])
rownames(merge4)<-as.character(merge4[,1])
dosage<-cbind(merge2[,5],merge4[rownames(merge2),5])
rownames(dosage)<-rownames(merge2)
merge2<-merge2[rownames(dosage)[which((dosage[,1]+dosage[,2])/2 >0.3)],]
merge4<-merge4[rownames(dosage)[which((dosage[,1]+dosage[,2])/2 >0.3)],]
